Polarization fields in nitride nanostructures: ten points to think about 
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Macroscopic polarization, both of intrinsic and piezoelectric nature, is unusually strong in III-V 
nitrides, and the built in electric fields in the layers of nitride-based nanostructures, stemming from 
polarization changes at heterointerfaces, have a major impact on the properties of single and mul- 
tiple quantum wells, high mobility transistors, and thin films. The concepts involved in the theory 
and applications of polarization in nitrides have encountered some resistance in the field. Here we 
discuss critically ten "propositions" aimed at clarifying the main controversial issues. 
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IIFV nitrides are a new frontier of semiconductor 
physics . One of their basic properties, macroscopic 
polarization, offers unique opportunities for device design 
and basic investigations. Polarization manifests itself as 
built-in electrostatic fields in polarized layers interfaced 
to each other. These fields affect the characteristics, per- 
formance, and response of multilayer nanostructured de- 
vices. Here we discuss some theoretical basics and a num- 
ber of practical implications of polarization fields for ni- 
tride nanostructures. 

The microscopic foundations of dielectric polarization 
theory, established in the first half of the present decade 
1^,^, identify the zero-field or transverse polarization Pt 
of a periodic bulk crystal as the gauge-invariant Berry's 
phase of its occupied Bloch orbitals, accumulated in an 
adiabatic transformation of the system from some refer- 
ence state to its actual state Pt, which can equiv- 
alently be viewed as the integrated polarization current 
flowing through the crystal during the transformation, 
has no relation with the charge density of the polarized 
dielectric. Importantly, Pt can be calculated accurately 
from first-principles density-functional calculations . 
The absolute polarization of a material can also be ob- 
tained by referencing its polarization to that of a system 
for which P=0 by symmetry (for wurtzite, this may be 
the zincblende phase Q). Although not directly mea- 
surable, Pt is relevant to polarization and fields in fi- 
nite, e.g. multilayer, systems through classical relations 
of the sort Pt = £^o Pl, where is the (measured 
or computed [^,0) static dielectric tensor, and Pl is 
the dipole moment per unit volume of the finite system, 
i.e. the termination-dependent longitudinal polarization 
(also identified with minus the screened field generated 
by the net polarization charge at the sample surfaces or 
outer interfaces). It is simple, as discussed below, to gen- 
eralize such expression to a multilayer systems relevant 
in technology in terms of the more fundamental Pt only 
In the rest of the paper, any mention of polarization 
refers to the transverse variant Pt- 



In the case of III-V nitrides 



we deal with the total 



polarization P — Pgp + Ppz(e) in a given strain state 
at zero temperature, in the absence of external fields. A 
key point is that in wurtzite nitrides the total polariza- 
tion at zero strain, known as spontaneous, is non-zero 
and large. The spontaneous polarization has a fixed di- 
rection on general grounds pO| , and specifically it points 
along (0001) in III-V nitrides. The strain-induced piezo- 
electric term Ppz(e) can in principle point in any direc- 
tion depending on strain. Since nitride multilayers are 
usually grown along the (0001) axis, Ppz(e) actually lies 
along that axis. Given the results of Ref. Q and the con- 
ventions explained therein, for a layer under tensile strain 
(lattice constant must expand to fit on substrate, e.g. Al- 
GaN on GaN) the piezoelectric polarization points along 
(0001), whereas for a layer under compressive strain (lat- 
tice constant must contract to fit on substrate, say InGaN 
on GaN) it points along the (0001) direction. Similarly 
to Pgp, the piezoelectric component Ppz(e) is also quite 
unusually large, given the giant piezoelectric constants of 
in-V nitrides H,|n|l. 



It is clear that the two components just discussed can 
sum up or cancel each other out depending on strain and 
polarity; whether this increases or decreases the field in 
the layer, depends on the polarization and geometry of 
the neighboring layers. In fact, the polarization of a ma- 
terial A manifests itself at the interface with a different 
medium _B as a fixed local charge accumulation. If the 
interface is insulating and gap-states-free (as is gener- 
ally the case for the isovalent common-anion nitride in- 
terfaces), simple relations holds |^,|l3| between interface 
charges and bulk transverse polarization. These were ver- 
ified directly in ab initio calculations on interfaces of bi- 
nary nitrides [^|^ and of other materials : polarization 
charge densities of order 10^'^ cm~^, localized in an inter- 
face region ~ 3 — 4 A thick are found to be common. As 
a consequence of these large charge accumulations [ p^ , 
macroscopic electrostatic fields exist in the interfaced lay- 
ers, screened by dielectric response [free-carrier screening 
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is taken up later on in the paper] . 

The typical effects of polarization fields are such, as 
discussed further below, that comparing the predictions 
of polarization theory with typical experiments is an in- 
herently indirect process, requiring considerable model- 
ing, ft is essential that this modeling be done appro- 
priately. We have previously applied a sophis- 
ticated (and complicated) approach to realistic nanos- 
tructures, whereby we solve selfconsistently an empirical 
tight-binding Schrodinger equation and Poisson's equa- 
tion. The latter reads 



dz dz 



dz 



Pi 



(1) 



where the position-dependent quantities D, e, V, p, n, 
and Pt are respectively the displacement field, dielectric 
constant, potential, hole density, electron density, and 
total transverse polarization. The effects of composition, 
polarization, and free carrier screening are thus included 
in full, and non-equilibrium as well as equilibrium carrier 
distributions can be studied. 

In this paper, instead, we discuss qualitative aspects 
of polarization effects referring to a simple expression for 
the the electrostatic field in the j-th layer of a periodic 
but otherwise arbitrary multiple quantum well or super- 
lattice made of layers of materials k of thicknesses Ik and 
dielectric constants Sk- 



Efc(Pfc-P,-)W£fc 



(2) 



with sums running on all layers (including the j-th). Fj 
is the total, spontaneous plus piezoelectric, transverse 
polarization of the material in layer j. This expression, 
which can be derived j|] under the somewhat restrictive 
assumptions of periodic boundary conditions and no free 
carrier screening, shows at a glance that polarization 
effects couple together all parts of a multilayer structure 
via polarization differences. Also, simply using Eq. |^ 
and the values calculated in Ref. for piezoconstants 
and intrinsic polarization, several often-neglected facts 
are immediately apparent. 

1 — Spontaneous and piezoelectric polarizations 
are comparable, and neglecting either of the two com- 
ponents leads to inaccuracies. This is best seen in nu- 
merical examples; for instance, Alo.2Gao.8N assumed to 
be strained on GaN has a piezoelectric component of 
-7.5 mC/m^ and a spontaneous one of -39.4 mC/m^; 
an Ino.1Gao.9N layer strained on GaN has -f 15.6 mC/m^ 
piezo, and -29.3 mC/m^ spontaneous. This reflects di- 
rectly on fields. For a strained GaN/Iuo.iGao.gN isolated 
QW (given the of GaN spontaneous polarization -29 
mC/m^) one gets a field of 30 kV/cm on neglecting the 
piezo component, while the full value is -1.63 MV/cm. 
For an unstrained isolated Alo.2Gao.8N/GaN QW, ne- 
glect of spontaneous polarization gives a field of -1.14 



MV/cm instead of the actual -1.89 MV/cm. The ef- 
fect is less dramatic in the AlGaN/GaN systems since 
the piezoconstants of AIN are extremely large so that 
the piezoelectric term is sizable even at low strain. Of 
course, besides the errors it causes, neglect of either of 
the components is conceptually ill-founded. 

2 Only polarization differences at interfaces are 
relevant. This fact is a consequence of electrostatics 
via V • P—-(J, it is consistent with the tenet of basic 
polarization theory that only polarization differences 
are meaningful, and is also borne out by Eq. ^, which is 
derived under simple electrostatics assumptions. Some 
workers incorrectly assumed earlier on, that in a generic 
MQW the field in the active (well) region was given by 
the polarization/dielectric constant ratio of the well ma- 
terial only. In the case of e.g. a GaN-matched MQW 
with equally thick layers of GaN and Alo.2Gao.8N, this 
assumption leads to predict a field of -3.2 MV/cm in the 
well, while using the correct procedure the fields are, as 
expected, equal and opposite in the well and barrier, and 
0.74 MV/cm in modulus. 

3 The geometry of the structure determines 
directly the fields. This comes about for two rea- 
sons. First, the field is determined by the structure 
geometry via Eq. ^; second, the field-induced potential 
drop changes with layer thickness. In a GaN-matched 
Alo.2Gao.8N/GaN system with 30 A wells and 60 A 
barriers the field in the wells is -1.31 MV/cm with a po- 
tential drop of 0.39 eV; if the wells are enlarged to 60 A, 
the well field is -0.98 MV/cm, and the potential drop is 
-0.59 eV. Hence, doubling the thickness of a single layer 
does not straightforwardly produce a doubled potential 
drop. This is most relevant since it is just this potential 
drop that causes the optical shifts visible in PL measure- 
ments. Note in passing that the number of barrier/well 
units and their growth sequence is also relevant. 

4 Piezo vs spontaneous dominance depends on 
composition. Nanostructures involving AlGaN/GaN 
will be dominated by spontaneous fields, due to the 
huge difference in spontaneous polarization between AIN 
and GaN, and their relatively small lattice mismatch. 
GaN/InGaN nanostructures will be dominated by piezo- 
electricity, since the spontaneous polarization difference 
between InN and GaN is small, whereas their lattice mis- 
match is large. This becomes obvious considering that, 
in the Vegard hypotesis, for each 10 % molar fraction 
change the specific strain is -1 % for InN in GaN, and 
+0.25 % for AIN in GaN, while AP^p is +0.3 mC/m^ 
and -5.2 mC/m^ for the same cases. 

The interplay of free carriers (intrinsic, extrinsic, in- 
jected, photogenerated, etc.) and polarization fields in 
nitride multilayers turns out to give rise to a rich variety 
of new phenomena. A closer look to these issues through 
accurate simulations and thoughtful experiments has in- 
deed revealed several hot spots: 



2 



5 Optical transitions are red-shifted, slowed 
down, and made less intense by the fields. This ef- 
fect, due to the progressive spatial separation of hole and 
electrons at the two sides of the QW, was observed by sev- 
eral workers in strained GaN/InGaN ||l^, and strained as 
well as unstrained AlGaN/GaN MQWs In calcu- 
lations, it was shown that at excitation densities typical 
for PL measurements, the shift in a single QW is perfectly 
linear with the thickness, and the oscillator strength is 
suppressed accordingly |jl^,|l3|,|l8| . The shifts are most 
often such that the transition occurs well below the band 
gap energy. Identifying the magnitude of the polariza- 
tion via measurement of the shifts requires an accurate 
modeling, expecially when different well thicknesses are 
used at the same time in a MQW 0, and also in the 
small-thickness limit where the field-induced red shift is 
not linear due to the competing confinement blue shift. 

6 Optical excitation, electrical injection, or dop- 
ing, lead to partial screening of the fields due to gen- 
eration of free carriers in the active layers or thereabouts. 
This was observed in experiment as a blue shift of the 
transition energy and recovering oscillator strength as a 
function of excitation power ||lj,|l^. As indicated theo- 
retically |l^,|l^, the blue shift is in fact relative to the 
field-red-shifted energy, and failure to recognize the red 
shift as due to polarization has lead to overestimates of 
the In contents in InGaN layers. Referring to Fig. |^, 
more points are worth a mention. 
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FIG. 1. Transition energies and probabilities of an isolated 
Alo.2Gao.8N/GaN quantum well as function of well thickness 
and excitation density (see Ref. Jist) 



A very substantial excitation is needed to recover ap- 
preciable oscillator strengths, and to blue-shift the transi- 
tion back towards the flat-band value; this finding agrees 
with thresholds of GaN-based lasers being up to one or- 
der of magnitude larger than in GaAs-based devices. Fur- 
ther, the oscillator strength never reaches unity, and the 
energy never equals the gap plus confinement energy even 
in the highest power limit. This is because the relevant 
length scale for free carrier screening of the field is of 
order 15-20 A due to the larger spatial extension of the 
quantum-confined carriers as compared to the polariza- 
tion charge. A possible exception are extremely thin wells 
where the screening length is similar to the well width, 
and the confinement blue shift competes with the field 
red shift. Note that while excitation-driven screening 
dies out as excitation is removed, permanent screening 
can be and has been achieved by remote doping ||l^,|l^ . 

7 The field-induced potential drop can become 
similar to the gap of the layer material, for typical 
fields (several MV/cm) over typical well thickness (~ 100 
A). Screening by the intrinsic carriers becomes impor- 
tant in this limit. As could be guessed, simulations show 
that at non-zero temperature, as a result of free carrier 
generation at the opposite ends of the sample, the field- 
induced potential drop across a given layer never exceeds 
the gap, but rather the field obeys |E| = Egi^p/d 
for thicknesses d > dc = Sgap/lEj. Since the relevant 
dc's are generally larger than typical QW thicknesses, 
the gap closure regime is generally not of interest there, 
but it may be to thick layers (used e.g. to measure alloy 
properties). On the other hand, cases exist where intrin- 
sic carriers determine the boundary conditions, as in the 
next point. 

8 — Self-doping may occur in "constrained" systems. 
This was observed experimentally by at least three 
groups [po|-p2| and confirmed theoretically |^3| in the in- 
stance of high mobility transistors. We pick the example 
of a thick (~ 300 A) AlGaN layer on a very thick (~ 
1.5 //m) GaN buffer, contacted at the top. The metal 
Fermi level, the Schottky contact barriers, and the inter- 
face band offset limit the maximum admissible potential 
drop through the AlGaN layer to about 1 eV. The po- 
larization field would instead imply a drop of about 4 
eV. For equilibrium to be established, a large part of the 
polarization charge at the heterointerface must therefore 
be screened by free carriers, coming in from the GaN 
layer (not from the contact due to the Schottky barrier) 
following the strongly attractive polarization-induced po- 
tential drop on the GaN side. The GaN layer is typically 
unintentionally n— doped to at least 10^^ cm^"^, hence a 
fair source of electrons. The end result at equilibrium 
p0|-p3| is the accumulation of a extremely high-density 
(up to 10^^ cm~^) two-dimensional electron gas at the 
AlGaN/GaN heterointerface channel (the mechanism is 
summarized, e.g., in Ref. p3| and discussed in detail 
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in Refs. |2gi, |2|l, and [g3|). The mobility and chan- 
nel transconductance are increased accordingly psf , also 
thanks to the absence of ionized impurities in the channel 
itself. Two groups [^,0 have interpreted their results 
in terms of piezoelectricity only, while a third p2[ has in- 
cluded spontaneous polarization. At the level of charge 
density CV profiling measurements, the uncertainties are 
sufficiently large, and the interpretations indirect enough 
that both ways of proceeding can get acceptable credit; 
as mentioned in the next point, however, only sponta- 
neous polarization stands the consistency check with the 
sample polarity. 

Two-dimensional hole gases can also be obtained 
thanks to fields, as demonstrated recently ||2^]. A su- 
perlattice is designed so that the field in the doped bar- 
riers and well is equal and opposite; this is achieved to 
a good approximation with a MQW with equal well and 
barrier thickness ( see Eq. ||). The field enhances, es- 
sentially by field ionization, the release of holes from the 
relatively deep acceptors in the GaN or alloy layer, and 
funnels the holes to the interface holding negative polar- 
ization charge. Holes accumulate at the interface forming 
a degenerate 2-D gas. Spontaneous polarization was ac- 
counted for, and measured charge densities matched well 
those predicted by modeling. 

9 — Polarity plays a key role: since the polarization 
in each layer has a definite sign dictated by strain and 
by the spontaneous polarization direction, the combina- 
tion of polarizations to produce fields (Eq. ^ is deter- 
mined by the polarity of the structure. This is in turn 
determined by the substrate- and preparation-dependent 
growth face [N-face, i.e. (0001) surface up, or rather 
Ga-face, i.e. (0001) surface up]. Hence control on, and 
awareness of the polarity issue is also essential, since it 
is easy to envisage how polarity influences experiments. 
The measurement of polarization effects through field- 
induced optical shifts in quantum wells is "sign-blind" , 
since the field sign is immaterial to the observed shift. In 
electrical devices such as HEMTs, if a charge accumula- 
tion is expected at a given interface for a given polarity, 
the occurrence of the opposite polarity will cause instead 
a depletion, as displayed in Fig. 0. 
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FIG. 2. Interplay of polarity and 2DEG formation in the 
HEMT channel (adapted from Ref. |22|). 

Of course, the magnitude of polarizations in the layers 
will determine the density of the 2-D gas, hence quan- 



titative agreement can be expected only including both 
components. Indeed, in CV experiments on controUed- 
polarity HEMT structures |^^, the appearance (or non- 
appearance) of a 2-D gas at a specific interface was shown 
to correlate |Q with the existence of the spontaneous po- 
larization component, and its value was estimated in fair 
agreement with our calculations Q|. 
10 — Piezoelectric and spontaneous polarization 
can be independently controlled to a large extent, 
basically because only the first depends on strain. It is 
easy to show [|l^ that it is possible to obtain completely 
unstrained, yet polarized layers, as well as strained yet 
unpolarized systems. In the first case, the piezo com- 
ponent vanishes by construction, hence spontaneous po- 
larization is unambiguously the only remaining source 
of field, and can therefore be measured. In the sec- 
ond case, the composition must be such that polariza- 
tion differences are zero across all interfaces, whence the 
fields vanish in all layers. It is rather clear that, since 
realizing this situation requires compensating piezo and 
spontaneous components, alloying with Al and In will 
be needed. For instance, the field-free case should be 
obtained with GaN-matched GaN/Alo.ioIiio.o6Gao.84N 
MQWs with equally thick wells and barriers. Also GaN- 
matched GaN/Alo.68lno.32N MQWs should exhibit the 
same behavior. In practice, one would measure the red 
shift vs. thickness for a series of GaN/AlInGaN MQWs 
starting at null In content; upon progressively increas- 
ing the In content, the shift should be first seen to van- 
ish, and then to reappear for larger In concentration. 
Problems (or opportunities) that should be watched out 
for in modeling are that the interface band offset may 
switch to type-II, and that alloy phase separation is likely. 
Of course, the "critical" concentrations mentioned above 
must be taken as indicative only, in view of their sensitiv- 
ity to the intricate interdependencies of strain, lattice pa- 
rameters, piezoconstants, and spontaneous polarization 
values. 
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